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Functions  Decreasing  in  Transposition  with  Applications  to 
Shock,  Damage,  and  Down  Time 

* » t 

M.  Hollander  , F.  Proschan  , and  J.  Sethuraman 

Florida  State  University,  Tallahassee,  Florida,  USA 

In  this  paper  a large  class  of  multivariate  densities  and 
functions,  including  the  multivariate  Poisson  distribution  and  the  com- 
pound multivariate  Poisson  distribution,  are  shown  to  have  the  decreasing 
in  transposition  property  introduced  by  Hollander,  Proschan,  and  Sethuraman 
(197T,  Ann.  Statist.  722-733).  Applications  relevant  to  reliability 
are  given,  including  applications  to  shock  models , cumulation  of  damage, 
and  component  down  times. 

1.  INTRODUCTION /ND  PRELIMINARIES 

In  this  paper  we  derive  a basic  theorem  (Section  2,  Theorem  2.1)  show- 
ing that  a large  class  of  mltivariatc  probability  densities  and  frequency 
functions  possess  the  decreasing  in  transposition  (DT)  property  (Defini- 
tion 1.3)  introduced  by  Hollander,  Proschan,  and  Sethuraman  [3]  (HPS[3])« 

Of  particular  interest  are  the  multivariate  Poisson  distribution  and  the 
compound  multivariate  Poisson  distribution  treated  in  Section  3*  Then, 
in  Section  4,  we  consider  applications  arising  in  reliability.  These  ap- 
plications pertain  to  shock  models,  cummulation  of  damage,  and  component 
down  times.  In  the  remainder  of  this  section  we  give  some  preliminaries, 
including  definitions  and  theorems,  which  will  be  useful  in  the  sequel. 

AMS  1970  subject  classifications:  Primary  62E99,  62H99;  Secondary  62N05 
Key  words  and  phrases:  Functions  decreasing  in  transposition,  multivariate 
distributions,  reliability,  shock  models. 
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HPS [ 3 ] define  a partial  ordering  on  n -dimensional  Euclidean  space 
(Rn)  presented  in  Definition  1.2  below. 

DEFINITION  1.1.  A vector  x = (x^,  . ..,  x^)  is  said  to  be  a simple 

transposition  of  a vector  x'  if  x and  x'  agree  in  all  but  two  coordinates, 

say  i and  j , i < j , x.  < x. , x!  = x. , and  x'  = x, • we  write  x 5 x' . 

ljij  J i*  — — 

Thus  x'  is  obtained  from  x by  performing  an  inversion  of  a single 
pair  of  coordinates  that  occur  in  their  natural  order  in  x. 

DEFINITION  1.2.  Let  x and  x'  be  two  n-dimensional  vectors  such  that 

there  exists  a finite  number  of  vectors  x°,  x1*  . . . , xk  in  Rn  satisfying 

x =.x°  £ ^ • ••  ^ xk  = xr;  i«e.»  x'  is  obtained  from  x by  a finite 

number  of  simple  transpositions.  We  say  that  x'  is  a transposition  of  x. 
HPS  [3]  then  define  a class  of  functions  as  follows: 

DEFINITION  1.3.  Let  A = (A, , ....  X ),  where  A.  £ A.  .£...£  A 

— 1 n 1 2 n 

are  n ordered  values  in  R1.  We  say  that  f(A^  x)  is  decreasing  in  trans- 
position (DT)  if: 

(a)  fCA*,  x")  = f(l_,  x)  for  each  permutation  it  = (ir  , ...»  tt^)  of  the 
indices  1,  2,  ...,  n,  where  A 11  = (A  , ...,  A ) and 

— TT_  TT 

1 n 

tt  , . 

x B ( x , . . . , X ) • 

n . 7T 

1 n 

(b)  x5x'  implies  that  fU,  x)  i f(A.t  X'). 


The  familiar  concepts  of  majorization  and  Schur  functions  will  be 
used  in  this  paper,  and  for  completeness  we  recall  their  definitions. 


DEFINITION  1.4.  Let  x^-j  2 ...  2 x^-j  be  a decreasing  rearrangement 
of  the  coordinates  of  the  vector  x.  Let  x and  satisfy: 


{ x[i]  * ^ y[i]»  i ~ n " 

i=l  i=l  L1J 


Jix[i] = J/ur 


Then  x is  said  to  majorize  £ (we  write  x ? jr ) . 

n 1 

DEFINITION  1.5.  A function  f from  R into  R is  said  to  be  Schur- 
convex  (Schur-concave)  if  x ^ implies  f(x)  2 (^)f(jr). 

HPS[3]  derive  basic  properties  of  DT  functions  and  show  their  rela- 
tionship to  other  classes  of  functions.  They  show  (HPS[3l,  Lemma 
2.2)  that  when  f(A^  x)  is  of  the  form  g(A^  - x ) , f(A.»  x)  is  a DT  function 
if  and  only  if  g(y)  is  a Schur-concave  function.  They  also  show  (HPS[3], 
Theorem  3.2)  that  the  DT  property  is  preserved  under  positive  mixtures 
and  (HPS[3],  Theorem  3.6)  that  the  DT  property  is  preserved  under  products 
of  positive  DT  functions.  They  also  prove  a "composition"  theorem  for 
DT  functions  (HPS[3],  Theorem  3.3)  and  a "preservation"  theorem  (HPS[3], 
Theorem  3.7)  for  Schur-concave  functions  under  an  integral  transform 


where  the  kernel  is  DT. 


. 
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By  restricting  ggty.*  z)  in  the  HPS  "composition”  theorem  to  be  of 
the  form  g(y_  - z)  , where  g(w)  is  Schur-concave , we  have  the  following. 


THEOREM  1.5.  Let  g^(x,  y_)  be  a DT  function  and  g(w)  be  a Schur- 
concave  function  such  that 

f (x,z)  = Jgj/x*  £)g(£-z)  dpfc) 

is  well-defined,  where  p is  a positive  permutation  invariant  measure. 
Then  ftx  , z)  is  a DT  function. 


Theorem  1.5  generalizes  the  following  result  of  Marshall  and  OH  in 
([5]*  Theorem  2.1). 


THEOREM  1.6.  ( Marshall -Olhin , [ 5 ] ) - Let  g^fx)  and  g2(x)  be  Schur- 

concave  functions  such  that 

f(x)  = Jg-^x  - £.)g2(z)dp(£) 

is  well-defined,  where  p is  as  in  Theorem  1.5.  Then  f(x)  is  a Schur-concave 
function . 

2.  DT  PROPERTY  OF  OVERLAPPING  SUMS 
In  this  section  we  state  and  prove  our  main  theorem.  Let,  for 
k — 2,  3)  ...,  n. 


= {X:  A is  a subset  of  size  k from  {1,  .... 


n}}. 


(2.1) 


k = 2 


* • • • * 


THEOREM  2.1.  Let  ” = O’  ....  X ) , {X,, 

n A 

be  independent  collections  of  random  variables, 
function.  Let  the  random  variables  in  {X^,  A e 


1 ' V> 


n» 


Let  X have  a DT  density 
A^J  be  i.i.d.  and  have 


a common  log-concave  density  function  g^,  k = 2,  ...,n.  Let,  for 
i * 1,  . . . i n, 

vs 

\ * x.  + V Jx*  . (2.2) 

k=2  A:  AcA^  and  i<A 


Then  Z ■ (Z^,  ....  Zq)  as  a DT  density  function. 


REMARK  2.2.  Mote  that  the  summands  appearing  in  the  Z, , ....  Z overlap 

l n 

considerably.  For  example,  appears  in  the  expressions  for  Z^  and  Z^, 

X122  appears  in  the  expressions  for  Z1,  Zg,  and  ly  etc.  Thus  the  inheri- 
tance of  the  DT  property  of  from  that  of  X is  complicated  by  the  over- 
lapping of  the  XV  . 

To  prove  the  main  result , ve  shall  find  it  helpful  to  have  available 
the  following  lemma: 

LEMMA  2.3  Let  k i.  2.  Let  {X^,  ’ c A^,}  be  i.i.d.  j andom  variables 

with  a common  log-concave  density  function  g (with  respect  to  the  counting 
measure  on  a lattice  or  - * '■  :r~  rvre).  Let 

- lxx  , i = 1 n.  (2.3) 

A:  AeA^  end  ieA 

Let  f(w,  , ....  w ) be  the  density  function  of  W . . ...  w ) . Then  f 
In  ~ 1 

is  a Schur-concave  function. 
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Proof.  We  will  now  give  a proof  for  the  case  where  g and  f are  den- 
sity functions  with  respect  to  counting  measures . The  proof  for  the  case 
where  g and  f are  density  functions  with  respect  to  Lebesgue  measures  is 
similar  and  will  be  omitted.  s 

Notice  first  that  ....  WQ  are  exchangeable  and  hence  f^,  ....  wq 


is  permutation  invariant. 


Fix  w 


w and  define 
n 


— /t J 


To  show  that  f is  Schur-concave , we  must  show  that 


P(A  ) s P(A  , ,) 

T.T  « * U'  .V 


VW2 


W1’W2 


(2.U) 


whenever 


i.e.  whenever 


(w1,  w 2)  ? (w^,  w£), 


0 £ B * a £ w 


(2.5) 


where 


w^=w  + a,  w2  = w -a,  w^=w  + B » w^  = w - 6 • 


The  following  desses  of  subsets  of  {lf  ...»  n),  will  help  us  to 
divide  the  elements  of  A into  those  that  contain  both  1 and  2,  and  those 
that  do  not  contain  either  1 or  2,  and  those  that  contain  exactly  one  of 
1 and  2.  Define 

V = {vs  v is  a subset  of  size  (k  — 2)  from  (3»  •••»  n))» 

U * {u:  u is  a subset  of  size  (k  — 1)  from  (3»  ...»  n}}* 


* 
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and 


S = {s:  s is  a subset  of  size  k from  {3,  ...»  n}}. 


Any  A in  is  of  the  form  12v,  lu,  2u,  or  s.  Fix  numbers  *12v,  v e V,  xg, 

s e S and  x , u e U and  define  the  event  B as  follows: 
u’ 


B * <X12v  * T £ V’  Xs  ‘ V s e S-  \u  * X2u  ' 2V  “ * UK  (2'6) 


Let 


2Y  = X.  - X_  , u e U. 
u lu  2u * 


Then  the  density  function  of  {Y^,  u e U}  conditional  on  B is 


vuIU 


* yu)fUu  - yu> 


y n[f(x  + y )f(x  - y 

Lnu  u u u u 


)] 


(2.7) 


and  the  conditional  probability  of  A given  B is 

l’w2 


P(Aw  w lB)« 
l’W2 


Inf(xu  + yu)f(xu  - yu) 


Tnf(x  + y )f(x  - y,  ) 

LU  U U 'u  u 


y ,u  "J 

'u* 


ifwj=  ?x12v+  K + Ix6*  J = 3,  ....  n 
v:  Jfv  u:  Jeu  s:  Jes 

and  w*  = Ix12y  + Ixu 


V u 

otherwise . 


Now,  for  each  xu,  f(xy  + yu)f(xu  - is  symmetric  and  log-concaye  in  y^. 

The  numerator  in  the  expression  for  P(A  |b)  is  therefore  a convolution 

1 ’2 
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of  symmetric  log-concave  densities  and  hence,  from  Theorem  1 of  Karlin  and 
Proschan  [V],  is  symmetric  and  log-concave  in  a.  Thus  if  (2.5)  is  satisfied, 
we  obtain 

P(A  IB)  S P(Av,  IB). 

1 ’2  1 * 2 


Inequality  (2.U)  now  follows  by  unconditioning,  and  thus  we  have  esta- 
blished that  f is  schur-concave.  | | 

We  now  give  the  proof  of  the  main  theorem. 


Proof  of  Theorem  2.1.  Let 

X^k)  = IXX  , i = 1,  ...  n,  k = 2,  ...»  n. 

X:  XeA^  and  ieX 

Let  X^  = (x|k\  ....  X^),  k * 2,  ....  n,  and  X = (X^  ...,  XR).  Then 

Z = X + X^2^  + + 

From  Lemma  2.3,  the  density  functions  of  ...»  X^n^  are  all  Schur- 

concave.  Since  the  density  function  of  Z is  the  convolution  of  a DT 
density  and  several  Schur-concave  densities,  it  follows  from  Theorems 
1.6  and  1.7  that  the  density  function  of  Z is  DT.  | | 

3.  APPLICATIONS  TO  MULTIVARIATE  POISSON  DISTRIBUTIONS 
In  this  section  we  present  several  applications  of  Theorem  2.1. 


3.1.  MULTIVARIATE  POISSON.  The  multivariate  Poisson  distribution 
(cf.  Teicher  [5],  Dwass  and  Teicher  [1])  can  be  defined  as  follows.  Let 
X^  be  a Poisson  random  variable  with  mean  nn  , i = 1,  ...,  n,  X^  be  a 


Poisson  random  variable  with  mean  m^. 


Ac , k—  2,  ..., 


n 


where  is 
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given  by  (2.1).  All  the  Poi-sson  random  variables  are  assumed  mutually 
independent.  Then  (Z^,  ...»  Z^),  as  defined  by  (2.2),  is  said  to  have 

a joint  multivariate  Poisson  distribution. 

Nov  we  specialize  further  to  the  case  where 

* 

m^  — m^ , Xc , k — 2,  ...,n*  (3.1) 

That  is,  all  X's  with  the  same  number  (>  1)  of  subscripts,  have  a common 
mean.  Then,  from  Theorem  2.1,  we  conclude  that  the  corresponding  frequency 
function  is  decreasing  in  transposition.  That  is,  let  m = (m^,  ....  mn) 

and  let  f(m,  z)  denote  the  joint  probability  mass  function  of  Z, , ...,  Z 

x n 

(where  the  other  means  are  held  fixed  as  in  (3.1)).  Then  when  m^  £ ...  s mn, 
and  z^  £ z/  , we  have  f (m,  z)  £ f (m,  z_'  ) . 

3.2.  COMPOUND  MULTIVARIATE  POISSON.  The  compound  multivariate 
Poisson  distribution  described  below  generalizes  the  compound  bivariate 
Poisson  distribution  considered  by  Holgate  [2]  in  the  context  of  certain 
ecological  situations.  Let  Z^,  Z^,  ...,  denote  the  number  of  individuals  of 

type  1,  2,  ...,  n,  respectively,  in  a quadrat  of  land.  V.’e  suppose  that  these 
individuals  arise  from  independent  clusters  and  assume  that  the  number  of  clus- 
ters N is  a Poisson  random  variable  with  parameter  X.  In  cluster  j,  there  are 

4 = 4 + i 14 

k=2  X:  XeA^  and  ieX 


individuals  of  type  i,  i = 1,  ...,  n. 
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where  X^,  X*|,  X e k = 2,  n,  are  independent  Poisson  random  vari- 
ables witn  parameters  m. , m. , k = 2,  n,  respectively,  nd  satisfying 

1 A 

* 

m^  = m^  for  X e A^,  k = 2,  n,  for  each  cluster  j.  We  have  already 

stated  that  the  clusters  are  assumed  to  be  independent . Thus 

= (Z^,  Z^)  has  a multivariate  Poisson  distribution.  The  vector 

Z_  of  the  number  of  individuals  of  the  different  types  in  the  quadrat  is 
given  by 

1 N 

Z = Z + ...  + z“ 


and  is  said  to  have  a compound  multivariate  Poisson  distribution.  We  now  show  that 
the  density  of  Z is  DT.  Conditional  on  N = N^,  Z is  the  sum  of  i.i.d. 

multivariate  Poisson  vectors  and  therefore  is  multivariate  Poisson  from 
the  infinite  divisibility  of  the  multivariate  Poisson  established  by.Dwass 
and  Teicher  [l].  Thus  the  conditional  density  of  !Z,  given  N = NQ,  is 


DT.  By  unconditioning,  the  density  of  Z is  also  DT.  This  proof  also 
shows  that  the  random  number  of  clusters  N,  could  be  any  random  variable 
taking  values  on  the  positive  integers. 


U.  APPLICATIONS  TO  RELIABILITY  PROBLEMS 


In  this  section  we  show  sample  applications  of  Theorem  2.1  in  some 
commonly  occurring  reliability  models . 


U.l.  NUMBER  OF  SHOCKS  EXPERIENCED  BY  COMPONENTS  OF  SYSTEM.  Consi- 
der a system  of  n components  experiencing  successive  shocks  of  various 
types  over  time,  each  shock  affecting  one  or  more  components.  More  speci- 


al 


■HMk 


MMWiw 
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fically,  in  a fixed  interval  of  time,  let  be  the  number  of  shocks  af- 
fecting component  i alone  (i  = 1,  2,  ...,n),  X^  the  number  of  shocks 

affecting  components  i and  J simultaneously  and  no  other  components 
(l  £ i < J £ n),  X^g  the  number  of  shocks  affecting  all  n com- 

ponents simultaneously.  Let  (Z^,  . . . , Z ) be  defined  as  in  (2.2).  Then 
denotes  the  total  number  of  shocks  experienced  by  component 
i,i=l,  ...,n.  Assume  that  (X^  ...,  X ) , (X^,  a e A^},  k = 2,  ...»  n 

satisfy  the  conditions  of  Theorem  2.1.  Then  we  conclude  that  the  Joint 

frequency  function  of  (Z,  , ...,  Z ) is  DT. 

l n 

h.2.  CUMULATION  OF  DAMAGE.  We  consider  a model  similar  to  that  of 
b.l.  However,  now  we  study  the  damage  accumulated  by  each  component  as 
a result  of  the  shocks  simultaneously  affecting  the  various  subsets  of 
components.  Let  D^.  be  the  damage  accumulated  by  component  i as  a result 

of  shocks  affecting  component  i alone  (i  = 1,  ...,  n),  D be  the  damage 

— J 

accumulated  by  components  i and  J as  a result  of  shocks  affecting  compo- 
nents i and  ,)  simultaneously,  but  no  other  components  (1  £ i < j £ n),  ...» 
and  P n be  the  damage  accumulated  by  each  of  the  components  as  a re- 


sult of  shocks  affecting  all  components  simultaneously.  Assume  that  the 
Joint  distribution  of  D * (D^,  ...,  Dq)  has  a DT  density  function,  that 

D^,  \ e has  a common  log  concave  density  function  f ^ , k = 2 , ...,  n. 
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Define : 

Ti  = Di  + l [Dx  , i = 1,  ....  n. 

k=2  A:  AeA^.ieA 

represents  the  total  damage  to  component  i accruing  from  all  shocks 

affecting  component  i (and  possibly  affecting  other  components  simultane- 
ously) , i = 1,  . . . , n. 

From  Theorem  2.1  we  conclude  that  the  joint  probability  density  of 

T..  , . . . , T is  DT. 

1 n 

4.3.  COMPONENT  DOWN  TIMES.  Suppose  that  an  n-component  system  is 
subject  to  repair  (replacement).  If  any  component  fails,  the  system  fails 
and  the  repair  (replacement)  time  for  that  component  is  recorded.  Like- 
wise, if  any  subset  of  components  fails  simultaneously,  the  system  fails, 
and  the  time  it  takes  to  repair  (replace)  each  of  the  components  in  that 
subset  is  recorded.  For  example,  the  operating  record  of  a two-component 
system  is  exhibited  in  Fie.  4.1  below. 


Component  1 


Component  2 


1-  Repair  - Repair  Repair  | 

8 1| T 1 | 10 1 

| Repair  — Repair  — • Repair  | 

| 5 1 | — 8 1 | 7 } | 10  1 


System 


Down 


- Down  — Down 


— Down 


Time 


4.1.  Operating  Record  of  a Two-Component  System 


f 


■w 
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Figure  4.1  indicates  that  a total  of  4 repair  (replacement)  periods  occurred 
during  [0,  t]:  1 period  was  due  to  the  failure  of  component  1 alone,  1 
period  was  due  to  the  failure  of  component  2 alone,  and  2 periods  were 
due  to  the  simultaneous  failure  of  components  1 and  2. 

Let  IL  = time  devoted  to  repair  of  component  i due  to  shocks  af- 
fecting component  i alone  (l  < i < n),  R^  = time  devoted  to  repair  of  each 

of  components  in  X due  to  shocks  affecting  components  in  X simultaneously 
and  no  other  components,  for  each  X e A^.,  k = 2,  ...»  n.  Assume  that  the 

Joint  distribution  of  R^,  ...,  Rfl  has  a DT  density  function,  that  has  a 

common  log  concave  density  function  f^  for  X e A^,  k = 2,  . ..,  n.  Assume 

further  that  R,  {R^ , X e } , k = 2 , ...,  n are  mutually  independent.  Define 

n 

Si  = Ri  + £ ER^  j i * 1>  ...,  n. 

k=2  X:  XeA^ieX 

represents  the  total  time  spent  repairing  (replacing)  component  i as 

a result  of  shocks  affecting  component  i,  i = 1,  ...»  n. 

From  Theorem  2.1,  we  conclude  that  the  Joint  density  function  of 
, . . . , Sn  is  DT . 
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